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Abstract 

As yet, there is no underlying fundamental theory for the transplanckian regime. 
There is a need to address the issue of how the observables in our present Universe are 
affected by processes that may have occurred at superplanckian energies (referred to 
as the transplanckian regime). Specifically, we focus on the impact the transplanckian 
regime has on two observables, namely: dark energy and the CMBR spectrum. We 
model the transplanckian regime by introducing a 1-parameter family of smooth non- 
linear dispersion relations which modify the frequencies at very short distances. A 
particular feature of the family of dispersion functions chosen is the production of 
ultralow frequencies at very high momenta k (for k > Mp). We name the range of 
the ultralow energy modes (of very short distances) that have frequencies equal or 
less than the current Hubble rate Hq as the tail modes. These modes are still frozen 
today due to the expansion of the Universe. We calculate their energy today and 
show that the tail provides a strong candidate for the dark energy of the Universe. 
During inflation, their energy is about 122-123 orders of magnitude smaller than the 
total energy, for any random value of the free parameter in the family of dispersion 
relations. For this family of dispersions, we present the exact solutions and show that: 
the CMBR spectrum is that of a (nearly) black body, and that the adiabatic vacuum 
is the only choice for the initial conditions. 
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1 Introduction 



There is still no fundamental physical theory of the very early universe which addresses 
issues that arise from the regime of transplanckian physics. One of these issues relates 
to the origin of the cosmological perturbation spectrum. In an expanding Universe, the 
physical momentum gets blue-shifted back in time, therefore the observed low values of the 
momentum today that contribute to the CMBR spectrum may have originated from values 
larger than the Planck mass Mp in the Early Universe. This is similar to the problems that 
arise in trying to explain the origin of Hawking radiation in Black Hole physics. In a series 
of papers [[I], |2|, [3], f|, §], it was demonstrated that the Hawking radiation remains unaffected 
by modifications of the ultra high energy regime, expressed through the modification of 
the usual linear dispersion relation at energies larger than a certain ultraviolet scale kc- 
Following a similar procedure, in the case of an expanding Friedmann-Lemaitre- Robertson- 
Walker (FLRW) spacetime, Martin and Brandenberger in Ref. (see also p|, |[ |9|, |T0|, |TT|| ) 
showed that standard predictions of inflation are indeed sensitive to trans-planckian physics: 
different dispersion relations lead to different results for the CMBR spectrum. 

It is the lack of a fundamental theory, valid at all energies, that makes the model build- 
ing of the transplanckian regime very interesting. The main issue is how much are the 
known observables affected by the unknown theory. The apparently ad hoc modification 
of the dispersion relation at high energies is contrained by the criterion that its low en- 
ergy predictions do no conflict the observables. Specifically, in this paper we address two 
questions: a) can the transplanckian regime contribute to the dark energy of the universe, 
and b) how sensitive is the CMBR spectrum to energies higher than the Planck scale Mp, 
where our current physics theory is known to break down. 

We choose a family of dispersion relations for the frequency of the wavefunctions that 
modifies the behaviour of the field at the ultrahigh energies of the transplanckian regime. 
The dispersion relation has the following features: it is smooth, nearly linear for energies less 
than the Planck scale, reaches a maximum, and attenuates to zero at ultrahigh momenta 
thereby producing ultralow frequencies at very short distances. We name the tail that 
part of the dispersion graph of very short distances that contains the range of ultralow 
frequencies less or equal to the current Hubble constant Hq. It follows that the tail modes 
are still currently frozen. We calculate the energy of the tail modes in order to address 
the former question (a) and show that although the tail does not contribute significantly 
to the CMBR spectrum, it has a dominant contribution to the dark energy of the universe 
[|T2f . The energy density of the tail modes is of the same order today as the matter energy 
density. 

The second question (b) is motivated by the problem that in most inflationary models 
the present large scale structure of the Universe is extrapolated from a regime of ultra- 
high energies (known as the transplanckian regime) originating from before the last 60 
e-foldings of the exponential expansion. In Refs. @, || the authors have demonstrated that 
the problem of calculating the spectrum of perturbations with a time-dependent dispersive 
frequency can be reduced to the familiar topic of particle creation on a time- dependent 



background ||13|| . We will use their observation in what follows. They also conjecture that 



the observed power spectrum can always be recovered only by using a smooth dispersion 
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relation, which ensures an adiabatic time-evolution of the modes. By taking the frequency 
dispersion relations to be the general class of Epstein functions [|14]], we check and lend 
strong support to their conjecture. We present the exact solutions to the mode equation 
for the scalar fieldj] with a "time-dependent mass", and the resulting CMBR spectrum 
below. We show that the major contribution to the CMBR spectrum comes from the long 
wavelength modes when they re-enter the horizon. The spectrum is nearly insensitive to 
the very short wavelength modes inside the Hubble horizon. 

The paper is organized as follows: in Section 2, we present the set-up and formalism 
of our analysis. The family of dispersion functions, exact solutions to the mode equations 
of motion and the resulting CMBR spectrum (from the Bogoliubov method) are reported 
in Section 3. In Section 4, we calculate the contribution of the tail modes to the dark 
energy of the universe today. In this work, we have neglected the backreaction effects of 
particle production. This assumption is fully justified from the calculation of the energy for 
the transplanckian modes, in Section 4. Due to the dispersed ultralow frequency of these 
modes, the energy contained in that transplanckian regime is very small (10 -122 ptotai) , thus 
the backreaction effect is reasonably negligible |15|, || . We present our conclusions in Section 
5. 



2 The Set-Up and Formalism 



Let us start with the generalized Friedmann-Lemaitre-Robertson- Walker (FLRW) line- 
element which, in the presence of scalar and tensor perturbations, takes the form O, [18| 

ds 2 = a 2 (ri) {-drj 1 + [<% + h(r], n)QSij 

+hi(rj, n)%- + h gw (rj, n)Q i:j dx % dx 3 
n 

where 77 is the conformal time and a(r]) the scale factor. The dimensionless quantity n is 
the comoving wavevector, related to the physical vector k by k = n/ a{rj) as usual. The 
functions h and hi represent the scalar sector of perturbations while h gw represents the 
gravitational waves. Q(x % ) and Qij(x l ) are the eigenfunction and eigentensor, respectively, 
of the Laplace operator on the flat spacelike hypersurfaces. For simplicity, we will take a 
scale factor airf) given by a power lawQ, a{rf) = Irjc/r/l 13 , where (3 > 1 and \r] c \ = /3/H(rj c ). 
The initial power spectrum of the perturbations can be computed once we solve the time- 
dependent equations in the scalar and tensor sector. The mode equations for both sectors 
reduce |]19| , |20| , ^TJ to a Klein-Gordon equation of the form 
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a 
a 







(2) 



where the prime denotes derivative with respect to conformal time. Therefore, studying 
perturbations in a FLRW background is equivalent to solving the mode equations for a scalar 

1 These functions are known for having exact solutions to second order differential equations in terms of 
hypergeometric functions. 

2 It has been argued in M that the analysis extends to other laws for the scale factor. 
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field n related (through Bardeen variables to the perturbation field in the expanding 



background. The above equation represents a linear dispersion relation for the frequency 

u, 

u 2 = k 2 = %. (3) 
cr 

The dispersion relation of Eq. (|3|) holds for values of momentum smaller than the Planck 
scale. There is no reason to believe that it remains linear at ultra-high energies larger than 
Mp. Yet, nonlinear dispersion relations are quite likely to occur from the Lagrangian of 
some effective theory obtained by the yet unknown fundamental theory. Nonlinear disper- 
sion relations, similar to the ones we consider in this work, are known to arise in effective 
theories of: nonlocal condensed matter or particle physics models arising from non-canonical 
kinetic terms [22, 23 1; from the dissipative behavior of a quantum system immersed into 



an environment after coarse-graining p^fl ; or from effective theories with phase transi- 



tions, time-dependent mass squared terms or effective potentials p5j |26| , g7]] . Perhaps, 
trans-planckian models motivated by superstring theory |28| , [29| or a two-stage inflationary 
model |3(| are plausible. In the latter case, one could easily envision for example a scenario 



with the first stage of inflation occurring at energy scales above the Planck mass Q followed 



by a nonthermal phase transition |3T| . The preheating [[JT], |32| from the nonthermal phase 
transition then leads to the second stage of inflation below Planck energies. In the former 
case, the motivation comes from the common belief that the superstring theory is the one 
that describes or at least is valid at energies of the transplanckian period. Taking this idea 
one step further, we incorporate the concept of superstring duality (which applies at trans- 
planckian regimes) in our analysis by choosing a particular family of dispersion relations 
that exhibits dual behavior^, i.e. appearance of ultra-low mode frequencies both at low and 
high moment a0 k. 

Despite the above comments and possible approaches, we should stress that any model- 
ing of Planck scale physics by analogy with the already familiar systems is pure speculation. 
We lack the fundamental theory that may naturally motivate or reproduce such dispersion 
behaviour. Nevertheless, it would be instructive to derive these dispersion relations from 
particle physics and string theory, as a step towards understanding the physical nature of 
the model. 

In what follows, we replace the linear relation u 2 (k) = k 2 = n 2 /a(r]) 2 with a nonlinear 
dispersion relation u(k) = F{k). The family of dispersion functions F(k) for our model is 
introduced in section 3. These functions have the following features: they are linear for 
low momenta up to the Plank scale kc = Mp, taken to be the cuttoff scale, but beyond 
the cutoff they smoothly turn down and asymptotically approach zero whereby producing 
ultra- low frequencies at very short distances. Therefore, in Eq. (0), n 2 should be replaced 



3 There is no reason why inflation must only occur below Planck energies. In principle, inflation at ultra 
high energies is equally possible. 

4 For example, when compactifying superstring theory in a torus topology, of large radius R and winding 
radius r, the frequency mode spectrum is dual in the sense that R and r are related as r = l/R. This 
means that each normal mode with a frequency n/R, where n is an integer, has its dual winding mode 
with decreasing energy that goes like 1/r = R [|25| . 

5 We would like to thank A. Riotto for pointing this out to us. 
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by: 



n 



' ff = a( V fF(k) 2 = a( V fF[n/a 



(4) 



We will also consider the general case of non-minimal or conformal coupling to gravity by 
keeping some arbitrary, unspecified, coupling constant £. Then, the equation for the scalar 
and tensor perturbations, that we need to solve, takes the form 



" i 



n 



-(1-60- 



= 



(5) 



For future reference, we define the generalised co moving frequency as 



tit 



TV 



eff 



1-60- 



(6) 



The dynamics of the scale factor is determined by the evolution of the background infla- 
ton field 0, with potential V(<f)), and the Friedmann equation. In conformal time, these 
equations are: 



a" 
a 



-4 + ^(/-aV W ) 



i" + 2-0' + a 2 
a 



dV(</>) 



0. 
0. 



(7) 
(8) 



Most of the contribution in the perturbation spectrum comes from long-wavelength 
modes, since at late times they are non relativistic and act like a classical homogeneous 
field with an amplitude /i given by: 



( 1 



V27ra 3 



dnn 



\fj, r 



1/2 



(9) 



These are produced at early stages of inflation, thus they are very sensitive to the initial 
conditions. The correct vacuum stateQ is given by the solution to Eq. (|5]) as clarified below. 
But if the stage of inflation is very long, and the Hubble parameter H is not changing 
considerably, then a Minkowski-like vacuum is a good first order approximation for the 
initial vacuum state (rj — > — oo), with: 



-inr) 



2n 



(10) 



However, if H was greater at the early stages of inflation, before the last 60 e-foldings, 
and/or the inflationary stage is short, then one must solve the wave equation and find the 
solution that minimizes the energy [33, |34, 27]. This is the correct vacuum state of the 



system. Otherwise, if one a priori chooses the Minkowski vacuum to be the initial vacuum 



6 In j^] the authors argue that there are two vacuum states. The argument extends to the criteria for 
choosing the right vacuum out of the two. Here we show there is only one true vacuum state which reduces 
to the Minkowski vacuum only at a certain limit. 
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state describing the system, the resulting value of (/a 2 ) is considerably underestimated as 
shown rigorously by Felder et al. [ |35|j . 

As already shown in Ref. |J, Eq. (H) represents particle production in a time-dependent 
background. We will follow the method of Bogoliubov transformation to determine the 
spectrum. The correct initial condition for the vacuum state is the solution to the equation 
that minimizes the energy. Hence, if the 'time-dependent' background goes asymptotically 
flat at late times, then in that limit the wavefunction should behave as a plane wave: 

pn ^n— oo ; . ^ n V _ (f f ) 



As it is well known on this scenario, in general at late times one has a squeezed state due 
to the curved background that mixes positive and negative frequencies. The evolution of 
the mode function \i n at late times fixes the Bogoliubov coefficients a n and /3 n , 

a _> an P -i^n ut v _|_ @ n P +i^ ut v (10) 

with the normalization condition: 

\a n \ 2 - \f3 n \ 2 = I . (13) 

In the above expressions, VL™ and Q™ 1 denote the asymptotic values of Q n when 77 — > =Foo. 
The spectrum of particles per mode is then calculated with the conventional Bogoliubov 
method j34j . The number of particles created n and their energy density p are calculated 
by the following expressions 

' f - ,21012 



(n) = J dnn 2 \(3 n \\ (14) 



2n 2 a 
1 



2tt 



kdk / u{k)duj \(3 k \ 2 . (15) 



If the resulting Bogoliubov coefficient /3& of the particles produced has a (nearly) thermal 
distribution, we can conclude that the CMBR in our problem is that of a (nearly) black 
body spectrum. We introduce the class of Epstein functions as the family of dispersion 
relations in Section 3, and derive the CMBR spectrum from the exact solutions to the 
evolution equations. 

Meanwhile, for the special features of our choice of dispersion relation, the modes at 
very high momenta but of ultra low frequencies u>(k) are frozen for as long as the Hubble 
expansion rate H of the Universe dominates over their frequency. We refer to that as 
the tail of the dispersion graph. In Fig. @, for the dispersed u 2 {k) vs. k, the tail 
corresponds to all the modes beyond the point kjj, where kjj is defined by the condition 
^{kn) — Hq, where H is the Hubble rate today. It then follows that the tail modes are 



still frozen at present. We calculate the total energy of the particles by using Eq. (IBP, as 
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well as the frozen energy of the tail. Thus the energy of the tail is a contribution to the 
dark energy of the universe: up to present it has the equation of state of a cosmological 
constant term. However, through the Friedmann equation, if is a decreasing function of 
time because until now it has been dominated by the energy density of matter and radiation. 
Therefore, whenever H drops below the frequency uj of an ultralow frequency mode, this 
mode becomes dynamic by picking a kinetic term and redshifts away very quickly. Hence, 
when the dominant contribution to the evolution equation for H comes from the tail energy, 
the behaviour of those modes with equations of motion coupled to the Friedmann equation 
becomes very complex. It is hard to calculate at which rate H drops in this situation. If 
eventually, H drops all the way to zero, all the modes in the tail must have decayed. Their 
equation of state, when H becomes zero, is that of radiation. The reason can be traced back 
at their origin in transplanckian regime. It is well known that scalar perturbations produced 
during inflation do not contribute to the total energy. Thus the origin of this modes is in 
the tensor perturbations. In their physical nature they correspond to gravitational waves 
of very short distance but ultralow energyQ We calculate their energy today in Section 4. 

We would like to elaborate on yet another possibility, which has not been mentioned 
before in the literature, that can give rise to a similar dispersion relation: very short or 
very large distance physics may have a curvature different from the FLRW element of 
Eq. ([I]), e.g. a different scale factor. This becomes clearer when recognizing the strong 
relation between the time-dependent dispersion relation and the curvature given by the 
time derivatives of the scale factor a(t]). The basic argument is the observation that the 
modulated frequency fl^ in the wave equation contains the contribution of these two terms, 
as given in Eq. (H). Therefore, while keeping the generalised frequency invariant, changing 
the first term in can be viewed or attributed to changes in the second term, such that: 

a" A" 

"n = a ^ nonlinear, k ~ U _ 6£) ~ = A ^Unear,k ~ (1 _ ®0 ~T 7 (16) 

where A is the new scale factor at very short distances. Even in the conformal case £ = 1/6, 
when the term a" /a drops out of Eq. (0), the time dependent frequency ^nonlinear k can 
mimic a term proportional to a" /a at X(rj)kc <C 1. Thus, any modulation of the dispersion 
relation is equivalent to a change in the behavior of the time-dependence of the background 
(a.k.a., the scale factor /curvature). In other words, we could have introduced a different 
curvature at very short (or large) distances from the start instead of a dispersed frequency^]. 

3 Exact Solution and the CMBR Spectrum 

We will consider the class of inflationary scenarios that through Eqs. (|7|) and ([8]) has a 
power law solution for the scale factor a(rj) in conformal time, a{rj) = |r/ c A/|^, with /3 > 1, 

7 We thank S. Carroll for pointing this out. 

8 We are using this equivalence in a sequencial paper |36) to demonstrate that a different large scale 
curvature of the universe is not possible as it conflicts with the observed CMBR data. Therefore, trying 
to reinterpret the SNla data in the light of a possible different curvature for the large scale regions of the 
universe may be ruled out. 
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Figure 1: Shown is our family of dispersion relations, for j3 = 1 and representatives values of ei (solid 
lines). We have also included the Unruh's dispersion relation (dashed line) and the linear one (dotted line) 
for comparison. 



and the following Epstein function []14] for the dispersion relation: 

<>• = ") = »' (if? + TT? + (TtM • (18) 

where x = (k/kc) 1 ^ = A\r)\, with A = (l/|r] c |)(n//c(7) 1 / /3 . This is the most general expres- 
sion for this family of functions. For our purposes, we will constrain some of the parameters 
of the Epstein family in order to satisfy the features required for the dispersion relation as 
follows. First, imposing the requirement of superstring duality, in order to have ultralow 
frequencies for very high momenta, we demand that the dispersion functions go asymptot- 
ically to zero. That fixes 

e 2 = 0. (19) 

On the other hand, the condition of a nearly linear dispersion relation for k < kc requires 
that 

| + | = 1- (20) 

Still we will have a whole family of functions parametrised by the constant ei, as can be 
seen in Fig. [I]. 

With the change of variables rj — > u = exp(A\r]\), the scalar wave equation (|5|) for the 
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mode /i n becomes: 



with: 



where: 



d 2 u + -d u + V(u) 
u 



V{u) 



+ 



£3 



60 a" 



u 2 (l + u) u(u+l) 



u 2 A 2 a 



(Jfccrl*!)' ( g 



(21) 
(22) 

(23) 



In the case of conformal coupling to gravity, £ = 1/6, Eq. ( pTj) is exactly solvable in 
terms of hypergeometric functions |14| . This is a well studied case in the context of particle 



creation in a curved background [13|. Even if we are not in the case of conformal coupling, 
the contribution 

«: = w+i) (24) 

a I] 2 

is going to be negligible at early times (77 — > —00); at late times, it can be absorbed in the 
dispersion relation Eq. fll8|) redefining the constants e,. 

As explained in Section 2, the correct initial condition is the vacuum state solution that 
minimizes the energy. In the case where 62 7^ 0, this vacuum state behaves as a plane wave 
in the asymptotic limit 77 — ► —00, with f^ m - ) — > y/e^n. However, when 62 = as in our case, 
the correct behavior of the mode function in the remote past is given by the solution of Eq. 
(HI) in the limit 77 — ► —00. The exact solution which matches this asymptotic behavior is 
then given by: 



/i ( m )( r7 ) = C l 



1 +u\ d 1 , ,1 



d-b,l + 2d, 



1 + w n 



■it 



where C m is a normalization constant, and 

d = id 



+ e 3 



(25) 

(26) 
(27) 



At late times the solution becomes a squeezed state by mixing of positive and negative 
frequencies: 

a n -in° ut ri , Pn 



2tt° n ut 



+ 



2fi°«* 



(2? 



with \/3 n \ 2 being the Bogoliubov coefficient equal to the particle creation number per mode 
n, and f2° u< ~ J~t[n. Using the linear transformation properties of hypergeometric functions 
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37 



, we find thatf] 



cosh ir(d + b) 



cosli7r(<i 

If d is a real number (eg > —1/4), then we obtain: 



b) 



(30) 



-2nb 



\Pn\ 



2 sinh 2ixb 



(31) 



It is clear from Eqs. ( |30"D and fl3l"D that the spectrum of created particles is nearly thermal 
to high accuracy^, 

|/3n| 2 ^e~ 4 ^. (32) 

Thus, we can immediately conclude that the CMBR spectrum is that of a (nearly) black 
body spectrum. That means the spectrum is (nearly) scale invariant, i.e., the spectral index 
is n s ~ 1. This is consistent with previous results obtained in the literature P, |8|, |9|, p~0f| , 
when using a smooth dispersion relation and the correct choice of the initial vacuum state, 
as discussed above. In Refs. and ||, dispersion relations, that were originally applied to 
black hole physics |], H , were used in the context of cosmology. New models of dispersion 
relations were proposed by the authors of Refs. M and |T0]. Our proposal for a 1-parameter 



class of models has a significantly different feature from the above, namely: the appearence 
of ultra-low frequency modes in the transplanckian regime. The implications of such a 
behavior for high momenta on the production of dark energy are discussed in the next 
section. 



4 Dark Energy from the "Tail" 

In Section 2 we defined the tail as the range of those modes in the frequency dispersion 
class (originating from the transplanckian regime), whose frequency is less or at most equal 
to the present Hubble rate, H (see Fig. 2). It then follows that they have not decayed 
and redshifted away but are still frozen today. Since H has been a decreasing function 
of time, many modes, even those in the ultralow frequency range, have become dynamic 
and redshifted away one by one, everytime the above condition is broken, i.e. when the 
expansion rate H dropped below their frequency. Clearly, the other modes have long 

9 In the most general case, where £2 in Eq. Jl8|), it is obtained jlij ]: 







cosh7r(d + 


b -a) 






cosh7r(d — 


b-a) 



(29) 



with a = \fe2. Also in this case the spectrum of the fluctuations is nearly thermal, with the parameter d 
controlling the deviation from thermality. 

10 We remind the reader that we have neglected the backreaction effects during the calculation, based on 
the result of a small particle number per mode, in the high momentum regime (fc 3> Mp) and a very small 
energy contained in these modes. Clearly, the particle number per mode being small is consistent with the 
result of the exponentially suppressed, near-thermal Bogoliubov coefficient. See Section 4 for the energy. 
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Figure 2: The range of modes in the tail, fejj < k < oo, defined by Eq. (|35|). iJo is the present value of 
the Hubble constant. 

decayed into radiation and the tail modes are the only modes still frozen. They contain 
vacuum energy of very short distance, hence of very low energy. The last mode in the tail 
would decay when and if H = 0. When the tail modes become dynamic by acquiring a 
kinetic term (when u(k) > H), they decay away as gravitational waves (explained in Section 
2). The tail starts from some value kn which must be found by solving the equation 

uj\k H ) = Hl. (33) 

The range of the modes defining the tail is then for kn < k < oo. Their time-dependent 
behaviour when they decay depends on the evolution of H and is complicated because they 
contribute to the expansion rate for H . Thus, their equations of motion are coupled to the 
Friedmann equation. 

However we can calculate their contribution to the dark energy today, when they are 
still frozen, thereby mimicking a cosmological constant. We calculate numerically (using 
'Mathematical the range of the modes in the tail from Eq. ( |3"3"D and use this value for the 
limit of integration in the tail energy given by Eq. (|i5|). Below we report these results for 
the case of a scale factor with (3=1 but other values of (3 were also considered numerically 
and they produce an even smaller dark energy due to the extra suppression in the integral 
coming from the Bogoliubov coefficient (3\. Eq. (^) is a messy transcendental equation 
but the solution to that equation is crucial to the dark energy since the value kn is the 
limit of the energy integral. That is why we solved Eq. (|33|) numerically and replaced it in 
Eq. ( |ToD for the energy, using different representative values of the parameter e\. 

The energy for the tail is given by: 

I r°° r 

(Ptaii) = ir- 2 \ kdk u(k)du \(3 k \ 2 , (34) 
2,11 Jk H J 
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while the expression for the total energy is: 



(ptotal) = J j U ( k ) duJ \Pk\ • (35) 

The numerical calculation of the tail energy produced the following result: for random 
different values of the free parameters, the dark energy of the tail is p ta u = 10 _122 /(ei), 
times less than the total energy during inflation, i.e. ■ £MiL = 10 _122 /(ei) at Planck time. 

Ptotal 

The prefactor /(ei), which depends weakly on the parameter of the dispersion family e±, 
is a small number between 1 to 9, which clearly can contribute at the most by 1 order of 
magnitude. 

This is an amazing result! It can readily be checked by plugging in the dispersion 
expression for Uk, Eq. ([H]), in the integral expression of Eq. ( [34]) for the tail energy, then 
using as the limit of integration the value kn found by the condition in Eq. ([33]). This 
result can be understood qualitatively by noticing that the behavior of the frequency for 
the "tail" modes is nearly an exponential decay (see Eq. (|TTD), and as such dominates over 
the other terms in the energy integrand of Eq. 



u 2 (k > k c ) « exp(-k/k c ) , (36) 

Hence, due to the decaying exponential, the main contribution to the energy integral in 
Eq. ([34]) comes from the highest value of this exponentially decaying frequency, which is 
the value of the integrand at the tail starting point, k^ ~ O(Mp), i.e., 

O-S^S* 10-122 ■ (37) 

Due to the physical requirement that the tail modes must have always been frozen, the tail 
starting frequency ui{k H ) is then proportional to the current value of Hubble rate H (Eq. 

We suspect this result is generic for any scenario that features ultralow frequencies which 
exponentially decay to zero at very high momenta for two reasons. First, all the modes 
with an ultralow frequency uj < H will be frozen and thus produce dark energy. Secondly, 
their contribution to the energy may be small because of the following. Due to this kind 
of dispersion in the high momentum regime, the phase space available for the ultra- low 
frequency modes with u{k > k c ) gets drastically reduced when compared to the phase 
space factor in the case of a non- dispersive transplanckian regime, controlling in this way 
these modes contribution to the energy density. The result for the tail energy also means 
that the tail energy dominates today's expansion of the universe. Thus, at present, we 
can not tell a priori the evolution of these modes and the time when they may become 
dynamic. Only the solution to the equation for the modes coupled (strongly at present) 
to the Friedmann equation would answer the question as to whether H will continue to 
decrease whith time. If that were the case, then these tail modes would also eventually 
become dynamic and decay. However, we calculated the equation of state for the limiting 
case when H = 0. In this case, all the modes in the tail are dynamic. The calculation 
of the energy density of the tail in the dynamic case from Eq. (|15D confirms that the tail 
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decays in the form of radiation, as expected since their physical nature is that of gravity 
waves of very short distances (but ultralow energy), originating from tensor perturbations 
during inflation. 

The opposite case is also a possible outcome to the coupled equations. It's possible that 
the frozen modes of the tail will prevent H from dropping further below, in which case 
these modes will never decay. We have not solved these coupled equation yet, therefore we 
are just speculating on the two possible outcomes of that calculation. The solution is left 
for future work. 

At present, these modes, originating from the transplanckian regime, are behaving as 
dark energy of the same magnitude as the current total energy in the universe. This idea 
is then a leap forward to this longstanding and challenging problem of dark energy, for 
at least two reasons: first, inspired by superstring duality, it is very plausible to speak 
of scenarios with ultralow frequencies and very high momenta. The tail modes, that are 
frozen at present, provide a good candidate for the dark energy as our calculations show. 
Secondly, although smooth dispersion functions that model the transplanckian regime do 
not affect the CMBR spectrum, this regime still leaves its imprints in the contribution to 
the energy of the universe. This is a rich and currently underexplored area to consider with 
respect to the cosmological constant mystery. 

5 Conclusions 

In this work we investigated two phenomenological aspects of transplanckian physics: the 
issue of dark energy production, and the sensitivity of the observed CMBR spectrum to the 
transplanckian regime. For this purpose, a family of dispersion relations is introduced that 
modulate the high frequencies of the inflationary perturbation modes at large values of the 
momenta k for the transplanckian regime. The smooth dispersion relations are chosen such 
that the frequency graph attenuates to zero at very high k, thereby producing ultralow 
frequencies corresponding to very short distances, but it is nearly linear for low values of k 
up to the cuttoff scale k c = M P . 

We present the exact solutions to the mode equations and calculate the spectrum 
through the method of Bogoliubov coefficients. The resulting CMBR spectrum is shown to 
be (nearly) that of a black body. This calculation lends strong support to the conjecture that 
smooth dispersion relations which ensure an adiabatic time-evolution of the modes produce 
a nearly scale invariant spectrum. Further, we elaborate on the issue of the initial conditions 
to which the spectrum is highly sensitive and show that there is no ambiguity in the correct 
choice of the initial vacuum state. The only initial vacuum is the adiabatic vacuum obtained 
by the solution to the mode equation. On the other hand, we showed that the assumption 
of neglecting the possible back reaction effects of the tail modes on the inflaton field is 
reasonable and is justified by the result of the tail energy calculation of Section 4. Also 
the Bogoliubov coefficient obtained, Eq. (p9|) , is exponentially suppressed, so back-reaction 
does not become significant. We would also like to stress that due to the dispersion class 
of functions chosen, defined in the whole range of momenta from zero to infinity, the total 
energy contribution of the modes produced is finite, without the need of applying any renor- 
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malization/subtraction scheme. In a sense, the regularization-renormalization procedure is 
encoded in the class of dispersion we postulate^. 

The most exciting result of this work is the generation of dark energy in the observed 
amount for the present universe 112] . This has its origin at the transplanckian regime, due 
to the presence of the dispersed tail modes with ultralow frequencies equal or less than 
the current Hubble constant. The evolution of these modes is given by their equation of 
motion, and it depends on the Hubble rate through the friction-like term 3H fi n . On the 
other hand, the evolution of the Hubble rate, given by the Friedmann equation, contains 
contribution from the energy of these modes. But currently, the Hubble constant dominates 
over their frequency in the mode equation of motion. It then follows that the tail modes, 
up to present, are still frozen and have been behaving like a cosmological constant term. 
Therefore their energy is dark energy^. 

We have calculated numerically the energy of the tail during the inflationary stage for 
different values of the dispersion parameter e. The calculation showed that the tail energy 
was 10~ 122 (times a prefactor /(ei) which weakly depends on e\ and, for random values 
of the parameter, takes values between 1 to 9) orders less than the total energy during 
inflation. This result is true for the whole class of dispersion relations. We chose different 
random values of the 1-parameter dispersion family ei, and the numerical calculation shows 
that e\ influences the energy at the most by less than an order of magnitude. We did not 
need to do any tuning of the parameters and used the Planck scale as the fundamental scale 
of the theory. Clearly, at present the tail energy dominates in the Friedmann equation, if 
the ratio of its energy to the total energy (as found by the calculation) was 10~ 122 /(ei) 
during inflation. The tail thus provides a strong candidate for explaining the dark energy 
of the universe (T^|. We suspect that the above result of producing such an extremely 
small number for the tail energy without any fine-tuning (and by using Mp as the only 
fundamental scale of the theory), is generic for any dispersion graph with a tail. The family 
of dispersion relations that feature a tail, corresponding to vacuum modes of very short 
distances, was motivated by superstring duality [^|, [29X 

However, in Section 2 we made the observation that introducing a dispersion relation 
is equivalent to introducing changes in the curvature of the universe, at very short or very 
large distances while keeping the generalized frequency Q n of Eq. (||) invariant. It is quite 
possible that the dispersion relation for the tail modes results from a different curvature of 
the universe at very short distances. This is an important link and we use it in a sequential 



paper [36] to demonstrate that the SNla data can not be reinterpreted away by changing 



the large scale curvature of the universe. Although it is counterintuitive, since large distance 
would correspond to low energy theories, we show in |36| that any changes in large scale 
curvature would disagree with the observed CMBR spectrum. 

It would be interesting to know what happens to the tail and the Hubble rate in the 
future. After all, a model is useful in so far as it can make future predictions. Although 
conceptually it is straightforward to find out the answer, given by the solution to the coupled 



"Because of these two results, we do not have the problems mentioned in Refs. ]T^, |l6| when discussing 
trans-planckian physics. 

12 For a different mechanism of generating a constant energy density, through the backreaction of cosmo- 
logical perturbations, that mimics a cosmological constant term, see Refs. [fM BM. 
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equations (|), (0) and (H), technically it appears difficult to predict the future evolution of 
the tail modes. The technical difficulty lies in the fact that at present, the tail equation of 
motion is strongly coupled to the Friedmann equation for H since the tail energy dominates. 
The Hubble rate would continue to decrease only z/these modes decay, but these modes can 
decay only when the Hubble rate decreases below their frequency. It may be possible that 
the frozen tail will sustain a constant Hubble rate which in turn will not allow the further 
decay of the tail modes. It is also possible that H will continue to decrease in which case 
the tail modes will become dynamic and redshift away in the future. It is only the solution 
to the mode equations coupled to the Friedmann equation, that will provide the answer on 
whether the Hubble constant and tail will decay in the future or remain at their current 
value. We do not have this solution yet, and the work is left for future investigation. In 
addition, the equation of state, w{t) = (p/p), is an observable that will provide a test to 
the model |30| , especially with the new data coming in the near future from the SNAP pl~ 
and SDSS [42] missions. 

However, we know that currently these tail modes are frozen vacuum modes of ultralow 
energy but very short distance, thus their energy behaves like a cosmological constant en- 
ergy. We also know they become dynamic and acquire a kinetic term only when the Hubble 
rate drops below the frequency. And, if they decay, the product is radiation of gravity waves 
at very short distances since their physical origin is from the tensor perturbations during 
inflation (it is well known that scalar perturbations do not contribute to the energy). The 
condition for the decay of the last mode in the tail is fulfilled when H has dropped to zero. 

Many of our results, e.g. the dispersion family and the exact solutions together with 
the Bogoliubov coefficients, could be applied to the Black Hole Physics. The issue of 
transplanckian physics was originally raised in the Black Hole context with respect to the 
sensitivity of the Hawking radiation to the blueshifted, superplanckian energy wavepackets. 
Following a phenomenological approach, a few dispersive models were introduced [0-fl in 
order to introduce a bound on the blueshifted energies and check the sensitivity of the black 
hole spectrum. We have introduced a new, different family of dispersive models, that also 
gives rise to a thermal spectrum. The analytical results of our class of dispersion models can 
be applied to the black hole physics and reproduce the same thermal Hawking spectrum. 
There are many subtleties involved due to the different symmetries of the two scenarios, 
but these issues are beyond the purpose and scope of this paper. It is left for future work. 
However, if the Hawking radiation for this class of dispersions is again thermal, it lends 
strong support to Unruh's conjecture that black hole radiation is insensitive to physics in 
the far ultraviolet (trans-planckian) regime, being predominantly an infrared effect. On the 
other hand, our class of models departs from the previous ones in the asymptotic behavior 
at very high momenta, with the presence of an infinite "tail" of ultralow frequency modes. 
The tail feature and energy results, applied to a black hole case, may raise interesting 
issues, in particular with respect to the black hole's information loss paradox. 
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